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Abstract. 

j^ 4. The rates of axion emission by nucleon-nucleon (NN) bremsstrahlung are reconsidered 

by taking into account the NN short range correlations. The analytical formulas for the 
neutron- neutron (nn), proton-proton (pp) and neutron-proton (np) processes with the 
inclusion of the full momentum dependence of an one- and two- pion exchange nuclear 
potentials, in the non-degenerate limit, are explicitly given. We find that the two- 
pion exchange (short range) effects can give a significant contribution to the emission 
rates, and are temperature dependent. Other short range nuclear effects like effective 
nucleon mass, polarization effects and use of correlated wave functions, are discused 
as well. The trend of all these nuclear effects is to diminish the corresponding axion 
emission rates. Further, we estimate that the values of the emission rates calculated 
with the inclusion of all these effects can differ from the corresponding ones derived 
with constant nuclear matrix elements by a factor of ~ 24. This leads to an uncertainty 
, factor of ~ 4.9 when extracting bounds of the axion parameters. 
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1. Introduction 

One of the most appealing solutions to the strong CP problem is the prediction of 
axions, pseudo Nambu-Goldstone bosons associated with the spontaneous breaking of 
the Peccei-Quinn symmetry p]-[l]. Light axions might be emitted from all types of 
stars and might influence their evolution. In the lack of any experimental evidence, it 
is important to derive bounds on the axion properties, i.e. on their mass (m ), on their 
couplings to matter (g ai , a = axion and i = photon, electron, nucleon, etc.) and on the 
axion decay constant (f a ). These quantities are related by the equations: 

ltfGeV 

f a ~ 6eV (1) 

m a 

9ai 7— (2) 

la 
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where N q is the number of quark flavours. These bounds can be obtained either from 
direct experiments or from astrophysical and cosmological considerations. For example, 
from the duration of the neutrino pulse from SN1987A one finds m a > 3 x 10~ 3 eV 
[5], while the requirement that axions do not overdose the universe leads to m a > 10 -6 
[6]. Bounds on f a (m a ) can also be extracted from the axion emission decay rates from 
newly born neutron stars during their cooling. The temperature and density in the core 
of such stars, just after their formation, can be T ~ 25 — 80 MeV and p ~ (6 — 8)po, 
where po ~ 2.5 x 10 l4 g/cm 3 is the ordinary nuclear matter density. In these extreme 
physical conditions only axions can excape from the core, and the dominant process for 
the axion emission is the NN bremsstrahlung: iViV — > NNa (N = neutron or proton). 
The decay rates of axion emission from these processes were previously calculated by 
different authors. Using an one pion exchange potential (OPEP) as NN interaction, 
Iwamoto calculated these emission rates for the degenerate (D) limit [7]. Using the 
nuclear matrix elements (NMEs) computed by Iwamoto, Turner calculated the emission 
rates for the non-degenerate (ND) limit [8] as well, while Brinkmann and Turner have 
performed these calculations for arbitrary degeneracy using a constant term for the 
momentum dependent NMEs [9]. This constant term was obtained by neglecting the 
pion mass as compared with the nucleon exchange momenta, in the expression of the spin 
summed squared matrix elements (the high momentum limit) [9]. They concluded that 
the ND regime is a better approximation for the axion emissivities than the D regime, in 
the physical conditions characteristic for a newly born neutron star. The main difficulty 
in these previous calculations was the numerical evaluation of the multidimensional 
integrals appearing in the emission rates, when the full dependence of the NMEs on 
the nucleon momenta is included. In our previous work [10] we developed a method 
of computing the axion emission rates by NN bremsstrahlung, with the inclusion of 
the full momentum contribution from a nuclear OPEP. We considered all the nn, pp 
and np bremsstrahlung processes in both the D and ND regimes, and found that the 
finite-momentum corrections to the emissivities are quantitatively significant for the 
ND regime and temperature-dependent. However, these calculations have neglected 
the influence of the dense nuclear medium. Ericson and Mathiot [11] have investigated 
the modifications on the axion emission rates induced by short-range nuclear effects 
like the exchange of two pions between the nucleons (which can be simulated by the 
exchange of a (heavier) p- meson between nucleons), by many-body polarization effects 
and by the consideration of a nucleon effective mass. They found that under the physical 
conditions relevant for the axion emission from supernova SN1987A, all these effects can 
give significant corrections as compared to the simple estimations in terms of the OPEP, 
and obtained a much smaller lower bound for the axion decay constant. In this paper we 
try to improve the calculation of the axion emission rates by including the corrections 
induced by p-meson exchange between nucleons, which mimick the two-pion exchange 
effects. For that, we use our method developed in [10], which allows the inclusion in 
the numerical computation of the multidimensional integrals, of the full dependence on 
momenta of the nuclear potetial. Further, we estimate the uncertainty induced by the 
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short-range nuclear effects when deriving of bounds of the axion parameters. 

The paper is organized as follows: in section 2 we describe our method of calculation 
and obtain the theoretical formulae for the axion emission rates. In section 3 we present 
our results and comment on them, while the last section is devoted to the conclusions. 



2. Method of calculations 



The axion emission rate by NN bremsstrahlung is given by Fermi's Golden Rule formula 
(see for instance [9]) 



^aNN = (27r) Z 



ui d3p 



(27r) 3 2^ 



(I " P " -E a x (SxZ\M\ 2 )5\P)F(f) (3) 



(2vr) 3 2£ a 



where, and E{ (i=l,4) are the nucleon momenta and energies, while p a and E a are the 
corresponding axion quantities; S is a symmetry factor taking into account the identity 
of the particles (S=l/4 for nn and pp channels and S=l for the np channel). 

F(f) = /z/g(i — fs){l — fi) is the product of Fermi-Dirac distribution functions of 

2 

the initial (1,2) and final (3,4) nucleons. In the non-relativistic limit E t ~ m + and 
using the non-dimensional quantities [9] yi — fii/T (/Xj are the chemical potentials of the 
nucleons and T the temperature) and U\ = p 2 /2mT, the expressions of the Fermi-Dirac 
functions read fa = (exp Ul ~ m +1) . The degenerate (D) limit satisfies y » 1, while in 
the non-degenerate (ND) limit y << — 1. For the spin summed squared matrix elements 
(NMEs) we use the following expressions [12J: 

S x S|M| 2 = S x ^g 2 aN m 2 ( -L) M NN (4) 

where g a N is the axion-nucleon coupling constant (we assume the same constant for the 
axion coupling to neutrons and protons). For the nn and pp the momentum-dependent 
factors M NN read [12J: 

m ( ' k ' 2 V i ( ^ V i (l-/3)|k| 2 |l| 2 

MnU = 2 + |l|2, 2 + 771772— 2N/I1|2 ■ 2^ (5) 

y|k| +m 2 1T J \\\\ +m*J (|k| +m;)(|l| + m%) 
while for the np process 

( Ikl 2 V ( 111 2 V lkl 2 lll 2 

M " P = L|2 ' J +2 \W^) -^-^ T^T- ' ' |2^ 2 , (6) 

\Jk| +mf r J +mj r J (|k| +m^.)(|l| + m*) 

where k = pi — P3 and 1 = pi — p4 are the nucleon direct and exchange transfer 
momenta, respectively. The last (exchange) terms in the above expressions arise from 
the interference of two different reaction amplitudes. They contain contributions from 
the scalar product (k • l) 2 , which has been estimated ([9], [12]) by replacing it with 
its phase-space averaged value, denoted by j3. From kinematical constraints j3 =0 in 
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the D regime, while it is ~ 1.0845 in the ND regime [32]. Since the most important 
contribution to the emissivity comes from the ND limit, we perform our calculations 
only for this regime. For calculation, it is worth to estimate the contributions of the 
momentum exchange terms as compared with the pion mass in the physical conditions 
of axion emission. A typical nucleon momentum is (3mArT) 1//2 . At temperatures T 
~ few MeV the momentum exchange terms are only a few times larger than the pion 
mass term, and thus the pion mass cannot be neglected in the denominator of the 
NMEs expressions ((5) and (6)). Also, at small distances, below 2 fm, one expects that 
two-pion exchange effects become important. Their influence on the bremsstrahlung 
processes can be estimated by the exchange of a p meson between the nucleons when 
the effective mass of this particle is taken to be m* w 600 MeV [H]. When these effects 
are included, a typical term in the NMEs expressions is modified as follows: 





Ikl 


2 


Ikl 


2 _|_ m 2 



with C p = 1.67 [H]. 

For practical calculations we follow the procedure of Brinkmann and Turner [9] to 
derive the ND limit, by performing the transformation to the center-of-mass system: 

Pi + P2 Pi - P2 

p+ = ^— ; p- = ^— ; 

P3c = P3 + P+! P4c = P4 + P+ 

=> Pi = P+ + p -; P2 = p+ - p ; ps = p+ + ps c ; 

P4 = P+ + P4c P4c = -P3c (8) 

and by defining the dimensionless quantities: 

U( =^L( i= l,4); M± =^; *.= J!L, (9) 

P+P P+P3c P-P3c nn x 

COS71 = | n r; cos^c = n -; cos^ = n r ; (10) 

|P+I|P-| |P+HP3c| |P-||P3c| 

Following our method described in refs. [10], [TS], we express the NMEs 
contributions as a sum of three terms: a constant term, corresponding to the high- 
momentum limit when the pion mass is neglected as compared with the momentum 
transfer terms in the NMEs formula, and two separated correction terms representing 
the one- and the two-pion (p) contributions, when the full momentum dependence of the 
NMEs is taken into account. This way, we can better compare our results with the ones 
coming from previous calculations, where only the high-momentum limit contribution 
is considered j9]. We did this by expressing the NMEs in terms of the scalar quantities 
a s = k 2 + l 2 = 2d(x + z/2) and a p = k 2 l 2 = d 2 [x(x + 2)sm 2 7 + ^ 2 /4], with d = 2mT, 
x = 2u3 C and z = u/T. In this way, one can arrange that the integral over the nuclear 
part to be reduced in a quite accurate approximation (within an error of ~ 2%) into an 



integral which is independent of angles, and the integration can then be performed for 
the general case of the full momentum dependence of the NMEs. 

Finally, we expressed these NMEs in the following compact forms: 

S x S|M| 2 = (J-}' [(3 - p) - \MlS - \M£\*\ (11) 

for the nn and pp processes. 

The corresponding one- and two-pion exchange corrections have the following 
expressions: 

ATT _ E>7T f<2 

\M n \ 2 = m 2 — (12) 

A similar procedure for the np process gives: 

S X S|M| 2 = 256m ^aN ^ 4 [(1 + 2p) _ | M ^|2 _ |M 2.| 2] (14) 

where: 

I Turn |2 2 L 'np L '^ ^ -"np ^np^^ /. _\ 



and 



K;i 2 = wn\ 2 (i6) 



The coefficients appearing in the above formula are polynomials depending on m, 
T and and on variables w_ and u 3c . The full expressions of the coefh cents A^, Bfj, 
C*j, DJp C, D and E can be found in [TP] , while the expressions of the coefficients A 2 ^, 
B 2 n l, CH, D 2 n l, Ell, F 2 l, G 2 n l and H% are given in the Appendix A. We also used the 
notation = cos^iCos^ c + sinjisiwycCos^, with <p the angle between the vectors p + 
and p . 

Further, we could express the axion emission rates in the ND limit in a form with 
a similar structure as NMEs: 

ND = ND^ l fnn fnn 

ann ta ™ yv) \ (3-/9)10 (3-/3)/ / 

where 

6^(0) = 2.68 x 10-y an e 2 ym 2 - 5 T^(f/m n y (18) 

is the expression calculated by Brinkmann and Turner [9]. Jo, C are double integrals 
over U- and u 3c , and are easily to be computed, while I 2 ^ is a five- dimensional integral 
which needs a separate treatment. 



I = / a/ (w_w 3c )(w_ - u 3c ) 2 e 2u -du^du 3c (19) 
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x du _ dU3c m 

V (2m_ + 2u 3c + mi f ) y ' 

PCO PU ic PIT PTT p2lT 

Inn= / / / / (n_ M3c )M^-^3c) 2 e- 2 "- 



JO JO JO JO 



«TWT«^ J iJ ^ lA# (21) 

For the np process we find a final expression analogous to Eq.(17): 



TTT T2lT 

np i 

3io ~ 3 1 



C = 0°)(( 1 + 2 / 3 )-^-^) (22) 



with: 



d?(0) = 10.72 x lO^^m^K) 4 (23) 



The correction integrals 2£ and have the following expressions: 



?rm 2 roo 



C = ^# / / V^*0(«- - «3c) 2 e- 2u - x 



o jo 

4(1 + 2/3) (m_ + m 3 c) + (7 + 2/3)mx 
(2u- + 2u 3c + m 1 ) 2 



du-du3 C (24) 



-^np — Inn (25) 



3. Results 



We performed the integrals of the Eqs. (19), (20) and (24) without dificulty. To 
integrate the 5th-dimensional integral (21) we built a special numerical routine based 
on multidimensional Gaussian quadratures, which is fast and higly accurate. In order to 
better understand the error of neglecting the contributions of the momentum dependence 
of the spin summed squared NMEs, we plotted in Figs. 1-4 the temperature dependence 
of the relative axion emissivities ((e(0) — e)/e(0)) for the ND regime: for the nn/pp 
processes (Figs. 1-2) and for the np process (Figs. 3-4). We considered two values for 
the p meson mass: i) the usual "free" mass m p = 770 MeV, and ii) an effective mass 
m* = 600 MeV, recommended when the meson exists in a dense nuclear medium |14j . 
The lower curves from figures represent the case when only one pion exchange (OPE) 
effects are taken into account in calculation, while the upper curves represent the axion 
emissivities when both the OPE and the two-pion exchange (TPE) effects are included. 
One observes that the correction due to the TPE effects is quite important especially 
at higher T, while the total correction due to the inclusion of both effects is larger at 
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lower T. Indeed, at T ~ 75K the relative correction to the e(0) due to OPE effects is ~ 
30 %, while this correction becomes ~ 57% when TPE corrections are included as well 
(see Fig. 2). The influence of the TPE effects diminishes at lower temperatures but the 
contribution of both of them to the emission rates calculated in the high-momentum 
limit increases to ~ 97% at T ~ 25 K. Also, one can see that the TPE corrections 
depend, as expected, on the effective p meson mass, being larger for the p meson mass 
modified by the dense nuclear medium (m* = 600 MeV). Another feature that we found 
is that these corrections are temperature dependent. 

Finally, one remarks that the consideration of the full momentum dependence of 
an OPEP and TPEP has the general tendence of diminishing the axion emission rates, 
and the total relative error in computing the axion emission rates is between 30% at 
T ~ 75K and 97% at T ~ 25K. This means that the axion emission rates might be 
overestimated up to a factor of 2, if one neglects the nuclear effects due to the full 
dependence on momenta of the nuclear potetial which includes both OPE and TPE 
effects. 

However, in a dense nuclear medium, other nuclear effects should also be taken 
into account in calculation, as some of them are more important than those described 
above. These effects induce additional uncertainties in calculating the axion emission 
rates. For example, nucleon effective masses are important, since the emission rates are 
proportional to mjf. The calculation of effective nucleon masses in a nuclear medium 
was subject of many works (see for example [16] and the references therein). In this 
sense, at nuclear densities p ~ 3po, the effective nucleon mass can be half of its rest 
value: m e ^ ~ 0.5m7v- Thus, the emission rates can be overestimated by a factor of 
~ 6. If one takes into account both uncertainities (full momentum dependence of an 
OPEP and TPEP and effective nucleon mass) the emission rates can be overestimated 
by a factor of ~ 12. At high densities, the short range effects need also the calculation 
of transition amplitudes with correlated wave functions, instead of plane waves (Born 
approximation). The modified wave functions read: 

C W = 0°(r)/(r); f(r) = y/l - j (q c r) (26) 

where jo is the Bessel function and the momentum q c ~ 800 MeV/c fixes the scale of 
correlations. However, the modification in the nuclear potential due to this correction is 
rather small because of cancellation between ir and p exchange effects. So, this correction 
is not expected to induce sizeable uncertainties to emission rates and axion parameters. 
Other sizeable effects come from polarization effects of the ir and p propagators from 
A-hole excitations. A reduction of emission rates of a factor ~ 2 is expected when 
both the correlation function and the polarization effects are taken into account [llj. 
Summing up all these effects, we estimate an overall reduction of the emission rates by 
a factor of ~ 24, as compared with the calculation when only the high-momentum limit 
of an OPEP is considered. Since the emission rate is proportional to the square of the 
axion coupling constant, the overestimation of the value of the emission rate induces an 
uncertainty of a factor of ~ 4.9 when deriving the axion parameters. 
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Figure 1. The axion relative emissivities from nn,pp bremsstrahlung processes, in the 
ND limit, as a function of temperature T. The p mass is taken m p = 770MeV. The 
lower curve represents the case when only OPE effects are included in the calculation, 
while the upper curve represents the case when both OPE and TPE effects are included. 




Figure 2. The same as in Figure 1, but for an effective mass m* = 600 MeV 
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Figure 3. The axion relative emissivities from np bremsstrahlung processes, in the ND 
limit, as a function of temperature T. The p mass is taken m p = 770MeV. The lower 
curve represents the case when only OPE effects are included in the calculation, while 
the upper curve represents the case when both OPE and TPE effects are included. 
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Figure 4. The same as in Figure 3, but for an effective mass m* = 600 MeV 
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4. Conclusions 

In this paper we calculated the axion emission rates by NN bremsstrahlung by taking 
into account the full dependence on momenta of a nuclear potential where both one- 
and two-pion exchange between nucleons were considered. The TPE were simulated 
by the exchange of a p meson between the nucleons. We obtained analytical formulas 
for the nn, pp and np processes, in the ND regime, and perfomed the integrals over 
the momenta with a good numerical precision. We found substantial reductions of the 
emissions rates in comparison with the case when in calculation was considered only the 
high-momentum limit of the spin summed squared NMEs. The corrections due to TPE 
(short-range) effects are larger at higher temperatures, while at lower temperatures the 
influence on the emission rates is reduced. However, the inclusion of both OPE and TPE 
effects brings a significant contribution to the emission rates, being as well temperature 
dependent. At T ~ 25 K, it brings a correction of about 97% to the values of axion 
emission rates calculated in the "crude" approximation of constant NMEs. Also, one 
sees that the TPE corrections depend, as expected, on the effective p meson mass, being 
larger for the p meson mass modified by the dense nuclear medium (m* = 600 MeV). 

We discussed as well other short range nuclear effects like effective nucleon mass, 
polarization effects and use of correlated wave functions. The trend of all these nuclear 
effects is to diminish the corresponding axion emission rates. Finally, we estimated the 
overall uncertainty in computing of the emission rates related to all these nuclear effects. 
We found that the values of the axion emission rates calculated with the inclusion of 
all these effects can differ from the coresponding ones calculated with constant NMEs, 
by a factor of ~ 24. This leads to an uncertainty by a factor of ~ 4.9 when extracting 
bounds for the axion parameters. 
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6. Appendix A 

We present here a part of the coefficients entering in the correction matrix elements for 
the cases which we presented in the section 2. 



AH = 8C> 3 T 3 («_ + u 3c f (C p - 2) + 4C> 2 T 2 ( u _ + u 3c ) 2 (m 2 ^ - 2m 2 ,) (Al) 
Bll = -A8C p m 3 T 3 («_ + u 3c f (u_u 3c )^ (C p - 2) + lQC p m 2 T 2 




(A2) 



= 64C p m 3 T 3 {u_ + u 3c ) (u„u 3c ) (C p - 2) + 



+16C p m 2 T 2 (u„u 3c ) (C p m 2 n - 2m 2 p ) 



11 

(A3) 



Dll = —QAC p m 3 T 3 (u„u 3c ) 1 (C p - 2) ( A A) 

Ell = 8m 3 T 3 (u_ + m 3c ) 3 + 4m 2 T 2 (u_ + w 3c ) 2 (m 2 + 2m 2 ,) + 

+2mT (it_ + w 3c ) m 2 (m 2 + 2m 2 ) + m 2 m^ (A 5) 

F 2 l = 32m 3 T 3 («_ + u 3c ) 2 + 16m 2 T 2 («_ + m 3c ) 2 (m_m 3c )5 + 

+16m 2 T 2 (u_ + u 3c ) (m_u 3c ) 2 ( m 2 + 2m 2 ) + 4mT (u.u 3c ) * m 2 (m 2 + 2m 2 ) ( A6) 

G 2 ; = 96m 3 T 3 (u_ + m 3c ) (m_m 3c ) + 16m 2 T 2 (m_m 3c ) (m 2 + 2mJ) (A7) 

^ = -64m 3 T 3 ( M _ M3c )^ (A8) 
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